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ABSTRACT 


A theoretical model for the heat and mass transfer in a 
double-effect solar still has been developed. An algorithm is also 
presented to determine the therlnal behaviour of the still. The 
model developed incorporates a detailed representation of the 
losses from the surface of the still and uses hourly and yearly 
meteorological data for Delhi to predict its performance. The 
special features of the solar still investigated comprise a thin 
film of water flowing continuously over the absorber plate 
constituting the first stage of the still and another stream of 
saline water flowing over the top serrated surface of the 
transparent plate. This exposes a larger water surface to 
insolation as compared to a simple basin type still of the same 
dimensions and thus constitutes the second stage of the still. 
Besides, the entire unit is kept at optimum tilt from the 
horizontal to receive maximum amount of solar radiation. The 
effects of various design, climatic and operational parameters on 
the performance of the still has also been studied. The assumption 
that the air is saturated at all points inside the still has been 
relaxed to bring the model as close as possible to the reality in 
practice. It is seen that the double-effect solar still is 35 to 
40 percent more efficient than the single-effect still in 
distillate production. Besides, the efficiency of the still can be 
enhanced by varying the operational and design parameters. 



CHAPTER - 1 

INTRODUCTION 

1 . 1 ENERGY: CRISIS AND RENEWABLE SOURCES 

It is a well known fact that energy represents the source of 
activities of human beings and nature. It was the discovery of new 
sources of energy like coal, petroleum products and 
hydroelectricity which provided a tremendous impetus to the 
industrialization and it was the quantitative availability of 
these energy resources which mainly contributed to the pace of 
development of the now industrialized nations. But .unfortunately , 
the pattern of consumption in the first half of this century has 
been, rather, wasteful as no attempt seems to have been made to 
conserve energy with a foresight into the future. 

The oil embargo of the early seventies shifted the focus to 
the renewable sources of energy. There is hardly any disagreement 
on the importance of harnessing renewable sources of energy, 
especially solar energy in the immediate future. At a time when 
world is becoming more and more aware of the problems of increased 
pollution and mounting cost .this clean source of energy deserves 
a serious attention in the world energy plans. 

1.2 SOLAR ENERGY 

Sun [2], the source of solar energy is a sphere of internally 

hot gaseous matter. The rate of energy generation, due to fusion 

12 21 

reactions in the sun is around 12 x 10 Q ( 1 Q= 10 JOULES). 
Out of this only a very small portion, around 5300 Q enters the 
Earth’s atmosphere. Even this small fraction is huge compared to 




2 


the total world energy consumption of 0.3 Q per annum. However, 

all this energy, that enters the atmosphere ,does not reach the 

surface of the earth. While at a point outside the atmosphere, 

2 

the energy intensity is 1353 W/m ,the average amount received on 

2 

the earth’s surface is just about 690 V/m . This value may vary 
with location, seasons (i.e. day of the year) and the time of the 
day . 

It is observed that the sites lying between the latitudes of 
30° North and South of the equator have sufficiently high solar 
intensities even during the winter months. India with its extreme 
latitudes of 8° and 32°North is fortunately very well placed for 
solar energy utilization. Actual solar radiation measurements in 
various ESCAPC Economic and Social Commission for Asia and the 
Pacific) countries have been used to produce charts of total 
solar radiation per unit area per unit time. The attached chart of 
Figure 1.1 shows that India receives the maximum solar radiation 
which for India works out to be 2100 KW-h per year per square 
meter area. 

However, the main problems being faced by the solar energy 
utilization are its low intensities requiring large collector 
areas ,low overall system efficiencies and the need of auxiliary 
equipment due to its intermittent availability. 

1.3 IMPORTANCE OF SOLAR STILL 

Water is a basic necessity of man along with food and air; 
the importance of supplying hygienic potable/ freshwater 
(containing less than or about 500 parts per million of salt) can 
hardly be over stressed. Man has been dependent on rivers, lakes 
and underground water reservoirs for fresh water requirements in 




Fig. 1.1 Average solar radiation ( in KWh per year 

per square meter area ) for ESCAP countries 
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domestic life, agriculture and industry. However, use of water 
from such resources is not always possible or desirable on account 
of the presence of large amount of salts and harmful organisms. 
Besides this there are several regions on the surface of this 
earth e.g the deserts, which have inhospitable climatic conditions 
and have only brackish water resources . In such places fresh 
water will have to be provided not only for domestic use, but also 
for agricultural needs. 

The only inexhaustible sources of water are the oceans. 
Besides, there are also several, hitherto unused, big lakes, 
inland seas and underground natural reservoirs containing salt/ 
brackish water. The chief drawback, obviously, is the very large 
salinity of such water. One ©f the very attractive schemes to 
tackle the problem of water shortage is the distillation of such 
water resulting in desalination. The conventional distillation 
processes are not only energy intensive but are also uneconomical 
for not too large demands of fresh water. However the developments 
in the use of solar energy have demonstrated that it is ideally 
suited for desalination, when the demand for fresh water is not 
too large. The rapid escalation in the cost of the fuels has made 
the solar alternative more attractive; in certain remote areas, 
this may be the only alternative. The least that can be said of 
solar distillation (distillation of saline water by use of solar 
energy) is that it is a viable option for providing hygienic 
potable water for a single house or a small community in most 
places of the earth. Further, the development of green-houses has 
resulted in minimizing the water requirements for carrying on 
agriculture on a small scale; thus solar desalination can support 
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small scale agriculture in regions having only brackish water 
sources . 

1.4 REVIEW OF PREVIOUS ¥ORK 

Solar distillation has been in practice for a long time. The 
earliest documented work is that of the Arab alchemists in 1551. 
In their historical review on desalination of water, NEBBIA and 
MENOZZI in 1966 129] mention the work of DELLA PORTA which he 
published in 1589. 

The great French chemist LAVOISIER in 1862 used large glass 
lenses, mounted on elaborate supporting structures, to concentrate 
solar energy on the contents of distillation flasks. The use of 
Silver or Aluminum coated glass reflectors to concentrate solar 
energy for distillation purposes has also been described by 
MOUCHOT [26J.Thus it appears that in 19th century scientists were 
familiar with harnessing solar energy for distillation not only by 
direct exposure to the sun but also by concentrating sun’s rays by 
means of mirrors and lenses. 

The conventional solar distillation apparatus (commonly known 
as solar still), was first designed and fabricated in 1872 near 
Las Salinas in Northern Chile by CARLOS WILSON. This was a large 
basin type solar still, meant for supplying fresh water to a 
nitrate mining community . 

No work on solar distillation seems to have been published 
after 1880’ s till the end of the First world war. With the renewal 
of interest in solar distillation, several types of devices 
have been described, e.g.roof type, tilted wick, inclined tray 
etc. Use of metal coated reflectors as solar concentrators for 
application in solar distillation has been described by ILAUSCH; 
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PASTEUR [30] in 1928 also used several concentrators to focus 
solar rays onto a copper boiler containing water. 

ABBOT in 19381271 used cylindrical parabolic reflectors 
(aluminuo coated surface) to focus solar energy onto evacuated 
tubes containing water. Be also devised a clock-work arrangement 
to track the motion of the sun. 

During the second world war, MARIA TELKES developed air 
inflated plastic stills for the U.S. Navy and Air Force for use 
in emergency life-rafts. 

The next stage in the development of solar stills was to 
improve the operating efficiencies of the various type of solar 
distillation devices. Forced air circulation was tried within 
stills to enhance the vapor condensation rate. Several 
investigators have attempted to make use of latent heat of 
vaporization released by the condensing vapor in either designing 
multiple-effect systems or for preheating the brine to increase 
the output of the stills. Several large scale distillation plants 
and integrated schemes for combining electric power generation and 
desalination of water have also been suggested as a way of 
improving the overall operating efficiency of the plant. 

The basin type solar still ( also called as the greenhouse 
type , roof type , simple or conventional type ) is in the most 
advanced stage of development. Several workers have investigated 
the effect of climatic, operational and design parameters on the 
performance of such a still. In attempts to improve the efficiency 
of distillation process, a large number of models for solar stills 
have also been suggested by various investigators. 
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1.5 ANALYTICAL APPROACHES 

There have been nunber of attempts to solve the heat and and 
mass transfer equations, applied to solar stills, analytically and 
predict their performance. A few important contributions in this 
direction are discussed below: 


1.5.1 SINGLE BASIN SOLAR STILL 

Cooper[251 and Frick [283 considered the periodicity of solar 
insolation and ambient temperature in the analysis of the solar 
still. However, the major drawback in their analysis has been that 
in the Fourier expansion of the periodic quantities only the 
fundamental harmonics have been retained. Also, they considered in 
their analyses an ideal solar still with many assumptions. Nayak 
et al.[12] presented the following analysis to predict the 
performance of a solar still with more general periodic behaviour 
of solar insolation and ambient temperature. A schematic sketch of 
the still analyzed by him is shown in Fig. 1.2 . The energy 
balance conditions at the top cover, saline water and the 
absorbing surface were be written as 


“g X <JTg/ dt = TgX 


. .( 1 ) 


M X dT / dt = To X H +q **q ~q ”q 

w w 2 s ^w ^cw ^rw ^ew 


.( 2 ) 




. . .(3) 
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Fig. 1.2 Single Basin Solar Still 
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q = Radiative heat transfer from water to glass cover 

i w 


= ^w’^ ^ - Tg) ...(6) 

q = Total heat loss from cover to ambient 

ct 


= h, X ( Tg - T^) 


Energy transferred from the absorbing surface to the saline water 
is written as 

Ok = *>3 »llx=o' 

while the loss to the insulation from the absorbing surface is 
given by 


Oi„s= - '^1 

where d^Cx,t) is the temperature distribution below the absorbing 
surface. The temperature distribution 9 Cx,t) in the insulation 
and ground is governed by the heat conduction equation 


p. C.09/at) = K. (aVdx^) j=1.2 (9) 

subject to the boundary conditions 

« (x=L.t) = 9Cx= L.t) (10) 

t 2 

- K, (39 /ax)l K = - K,(a»,/ax)i K <ii) 

1 1 •X=L Z 2 'X=L 

and 


9 (x,t) is finite for x > ® (12) 

2 

The solar intensity and ambient air temperature were considered to 

be periodic in time and Fourier analyzed in the form 

H = a + Re y A„ exp { i(n(i)t - a ) 1 . . . . ( 13) 

so " n n 

n 

and 

T = b + Re y B exp [ i (n«t )1 ....(14) 

a 0 “ n n 

n 

where o = 2ir/(24x 60 x 60) sec * 

The following periodic solutions were assumed. 

«^(x,t)=A^ X +B^ +Re J] exp(-p^^ x) 

n 


for 0^ X :s; L 
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9 (x,t) = B + Re r C exp <in»t - 6 x) for x i L 

2 2 “ 2n 

n 

Tg (t) = + Re J] exp (in«t) 

and 

(t) = + Re j; exp (in«t) 

n 

where 

“t <1 

B = o, (1 + i) 

^2n 2 

a. = (a p. C,/ 2K j = 1 , 2 

Jf *1 «J 

The constants , B . A , B . C , g , g , H and H were 

1122 1 n 0 no n 

deternined from equations given above and the exact solutions 
obtained . 

1.5.2 SIMPLE MULTIPLE WICK SOLAR STILL 

Kumar and Sodha et al.f41 presented a steady state analysis 
of simple multiple wick solar stilK Fig. 1.3). Using Dunkle's 
relations and neglecting the heat capacity of the irradiated mass 
of water and the glass cover, the energy balance for the glass 
cover and the water sheet were expressed as 



M 

X 

+q 

s ^rw 

+q -►q 

^cw ^ew 

-■■a = » 

• • 

.(15) 

and 

To X H 

dt s 

"Q 

^rw 

-q -q 

^cw ^ew 

-“b X (Tj, 


. . . (16) 

where 

the terms 

q 

cw 

* ‘^ew ’ 

’rw 

q the 

& 

same as already 


explained in equations (4), (5), (6) and (7), respectively and 


l/hjj= L/K + 1/hj 

The amount of water distilled per unit time per unit area is 
obtained by using Dunkle’s relations. 

= Qe^ / X = (16.273 x lO’^x (P^ - Pg)l)/X (T^ - Tg)..(17) 

where the temperature of water T^ and that of the glass cover T^ 
were evaluated by the simultaneous solutions of the equations (IS) 
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and (16) by numerical methods. 

1.5.3 DOUBLE BASIN SOLAR STILL 

Tewari et al.[14] have analyzed the performance of a 
double basin solar still, usin^ the periodic theory developed 
earlier. The major heat flux components are shown in Fig. 1.4. 
The heat balance equations are 


‘gu 

dT^ / 

II 

H 

s 

+ 

% 

• ’a 

[ 

wu 

dT / 

wu 

dt = T 

2 

H 

s 


h 

5 

^ ■ ‘’u 

‘gL 


dt = T 

3 

H 

s 

+ 

^L 

- ■'a^ ’■gL - Twu^ 

'wL 


dt = 

Hs 

+ 

h 

3 

' »ix=a - \l? - Ol 

H 

5 

» ' ““a 

o 

II 

H 

- T 

wL 

) - 

K,<a9/ax)|^,„ 


where 

q^j = Total heat loss from water surface in the upper basin to the 

upper glass surface 

= q + q + q 
rwu cwu ewu 

qj^ = Total heat loss from water surface in the lower basin to the 
lower glass surface 

**rwL * *^cwL * **ewL 

u and L refer to upper and lower glass covers respectively and 
q I q 1 q and q , . q , . Q , are the same as a . a < 
q^^ explained earlier provided T^ and T^ are replaced by T^^ and 
Tg^ and T^j^ and T^j^ respectively. 

Energy balance at the surface of insulation in contact with air is 
given by 

- K^(99/ac)|^^j^= f »li = L - 

The solar insolation and ambient air temperature can be considered 
periodic and are given by equations (13) and (14) respectively. 





14 


The following periodic solutions were assuaed. 

«(x,t) =A X +B +Re 1 {C^ exp(-^^ x) +D^ exp(p^ x)) exp (inwt) 

^rfu “ ^0 ^ (inwt) 

TgL ~ ^0 * ^ (inwt) 

^wu ■ ^0 ^ **n Cinwt) 

n 

^wL “ ®o ■*■ ^ ®n (in»t) 

n 

where the constants A. B, C , D„, g^. g^, h . h , G . G . H and 

nnononono 

were determined from the above equations and the exact 
solutions obtained. 

1.6 PRESENT WORK 

In the works reported so far, it has been assumed that the 
air is saturated at all points inside the still. Besides, the 
design of the stills have almost been identical in the sense that 
all of them have been basin type stills in which a large volume of 
water is kept on the top of the glass plates resulting in large 
attenuation of solar insolation passing through the thickness of 
water layer and reaching the absorber plates 

The present analysis of heat and mass transfer phenomena in a 
double-effect solar still discards a very important assumption, 
that the air is saturated at all points inside the still, used in 
previous works in this field. An improved equation for 
transmission of solar radiation as suggested by P.T.Tsi 1 ingr is 
[20] has also been used. Besides, the design of the still being 
investigated is such that it exposes a larger water surface to the 
solar radiation as compared to the solar stills of same dimensions 
analyzed in previous works. There is a continuous flow of water in 
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the lower stage of the still over the absorber plate which is in 
contrast to the solar stills analyzed in earlier investigations. 

In the present work, a computer model of a double-effect 
solar still has been developed. The input to the program are the 
hourly and monthly solar radiation and ambient temperature data 
for New Delhi chosen as a representative North Indian location. 
Also information about the dimensions of the still, operating 
conditions and thermophysical properties of water and air are 
provided . 


CHAPTER “ 2 

MATHEMATICAL MODEL 


The double-effect solar st i 1 1 , considered in the present work, 
consists of three distinct plates: the top glass cover, the 
intermediate glass plate and the absorber plate. The thermal 
performance of such a still depends on the heat capacity of these 
three plates. 

Figure 2.1 shows a schematic diagram of the aforesaid solar 
still to be analyzed. The still consists of one glass cover at the 
top and one absorbing plate at the bottom, as well as one glass 
plate in the middlejall of these plates are parallel. In order to 
collect more insolation, the equipment is tilted from the 
horizontal, with the cover inclined towards the sun. Accordingly, 
several weirs are installed in parallel on the upper surface of 
the glass plate to hold the saline water, while a blackened wet 
jute cloth is placed on the upper surface of the absorbing plate 
to form the liquid surface. Except for the glass cover, all parts 
of the enclosure are carefully insulated with a thick layer of 
asbestos to make the heat loss as low as possible. During 
operation, the top glass cover and intermediate glass plate 
transmit solar radiation and the absorbing plate is then heated 
directly by solar radiation. The saline water feeds are introduced 
on the upper surfaces of both the absorbing and intermediate glass 
plates, where some water evaporates, while the remainder is 
collected at the bottom and discarded as concentrated brine. The 
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vapor produced from the saline water rises from the absorbing 
plate and is condensed on the lower side of the intermediate glass 
plate. The heat of condensation given up by the condensing vapor 
on the lower side of the glass plate is conducted through this 
plate and furnishes heat to evaporate an equivalent amount of 
water from the saline water sliding down the upper side of the 
plate. Finally, the heat of condensation given up by the 
condensing vapor on the lower side of the glass cover is 
transmitted through this cover to the ambient. The condensates 
produced from the first and second effects are collected as fresh 
water in a trough from the enclosure. 

2.1 DEVELOPMENT OF THE MODEL 

The still is considered to be one dimensional with respect to 
the heat transfer, with the temperature drop through the cover and 
the plate thickness being neglected. Besides, the air at the water 
surface may be considered to be saturated while that at the 
condensing surface may be unsaturated or supersaturated. This 
assumption made in the present study is a significant departure 
from the earlier theoretical analyses which have assumed that 
water is saturated at all points inside the still. The model 
considers the transmission of the incident solar radiation through 
the surface, its attenuation as it passes through the subsequent 
water layers and the final absorption in the absorber plate. The 
blackened jute cloth is assumed to be perfectly absorbing and the 
radiation incident on it is completely absorbed there. It is 
further assumed that there are no losses from the bottom and sides 
of the still, these being perfectly insulated from the 
surroundings. The temperature of the absorber plate being higher 
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than that of other plates, it looses heat to the interuediate 
glass plate which looses heat to the glass cover and the glass 
cover loosing heat to the surroundings by convection and 
radiat ion. 

The above mentioned processes are discussed in detail in the 
ensuing sections and presented in a mathematical fora, using basic 
principles of heat transfer. 

2.2 ATTENUATION OF SOLAR RADIATION IN A WATER LAYER 

The performance of a solar still largely depends on the 
amount of the radiation which reaches the absorber plate at the 
bottom of the still. 

The most widely used mathematical expression on transmission 

of solar radiation through a water layer in literature is the Rabl 

and Nielson 111] four term exponential expression given by 

4 

f(z) = li X I J exp (-n. z) ...CD' 

“ i = l * 

where and n^ are empirical coefficients and have different 

values for different wavelength bands. 

For water thicknesses in the range from 1 cm to 10.0 m a 
relation for the fraction of incident solar energy f(z) 
transmitted through a thickness of z cm was obtained by BRYANT 
and COLBECKIlOl ,as 

f(z)=a' - b' ln( z ) . . .(2) 

where a' =0.73 and b' =0.08 (a' =0.36 if z is in meters). This 

representation is simpler for use than that of RABL and NEILSON. 
However, for a water thickness <1.0 cm, Eq (1) is not accurate, 
especially for the region very close to water surface. For this 
reason, a modified expression [19] for f(z) in this particular 
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region is introduced as follows: 

f(z) = a - b ln( z+c ) ...(3) 

Here the constants at b and c of Table 2.2 have been determined 
from the experimental data given in Table 2.1. The comparison of 
Eq(3) with Eq(2).as well as with RABL and NEILSON’s expression, 
are presented in Figure 2.2. It is seen from this figure that 
Eq(2) is the most accurate for the range of water thickness from 0 
to 1.0 cm. 

2.3 ABSORPTIVITY AND TRANSMITTANCE OF A WATER LAYER ON A GLASS 
PLATE 

If the reflection of solar insolation at the surface of a 
water layer on the upper surface of a glass plate is neglected , 
then 

As shown in Figures 2.1 and 2.3, we consider the optimal design 
case in which the height of weirs, as well as the distances 
between the weirs, are so installed that the depth of water layer 
at one end is exactly zero. Thus, the depth of water layer is 
related to the y-axis by 

z = p X (1 - y/L) ...(5) 

Since the solar energy flux transmitted through the glass cover is 

(t_x I ), that transmitted through the water layer on the glass 
g 0 

plate may be calculated from 

(Xg xI^)xT^=(l/WxL)I(XgXl^) jf(z)x W dyj 

Bence, t = (1/WxL) I f(z)x W dy) (6) 

* Jo 

From Eq (5) we have 


dz = -(p/L) X dy 


(7) 
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Fig, 2. 2 Fraction of incident energy remaining after 
passing through a thickness z 
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Using the expression for f(z) from Eq(3) and using Eq<6), the 
expression in Eq(6) reduces to 

P 

= (1/WxL) / {a-b ln(z+c)) x (L/p) dz . . . . C8) 

o 

Performing the integration in Eq(8), we get 

P 

= (1/p) X {axp - bxp ln(p+c) + b / [z/(z+c)l dz } 

0 

= (1/p) X {axp - bxp ln(p+c) + bxp - bxc ln(p+c) + bxc ln(c) ) 

= <a-bx{ I (p+c)/pjxln(p+c)-(c/p)xln(c)-l }> ....(9) 

The solar energy flux absorbed by the water layer is obtained by 
using £q(4) and the fact that 

f(0) = 1 = a-bxln(c) 

Hence , 

“w ' “ - 

= b{ [ (p+c)/pl { ln(p+c)/c J-1) ....(10) 

2.4 (KIVERNING EQUATIONS 

2-4.1 ENERGY BALANCE FOR THE TOP GLASS COVER 

The energy exchange of the glass cover with the ambient and 

the transparent plate is shown in Fig. 2.4. The glass cover is 

assumed to be at an uniform temperature T . The glass cover 

o 

receives energy from the sun part of which is absorbed by it 
and some part of it is reflected while most of it is transmitted 
through it. 

Incident radiation energy absorbed = a^x 
Incident radiation energy transmitted through glass cover= 
Evaporation heat gain from intermediate glass plate= 

Convective heat gain from intermediate glass plate®Q^^ 

Radiative heat gain from intermediate glass plate=Q ^2 
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Convective energy loss from cover to ambient=Q =h (T -T ) 

'ca a c a 

where . the wind convective heat transfer coeff icientCT J , is 
given by 

h^= 5.7 + 3.8 X V , W/m^-K 

where V is wind velocity in m/s 

Energy lost by radiation to the surroundings=Q =e o(T ^-T 

ra g c a 

Using the above expressions, we can write the total inward flux of 
the glass cover as 


+ Q + Q + Q 

g 0 62 ^C2 ^r2 

and the total outward flux of the glass cover as 

Q + Q 
ca ^ra 

considering steady state performance of the still. 


o -X I + Q +Q +Q =Q +Q 
g o 62 C 2 ^r2 ca ^ra 


where 


we can write 
( 11 ) 


^2 ■< 

O'wgo’' V hweg « '^2- 

0.9x o X 

A^= Ajj/sin* 

A = A. xtan<b 
we b ^ 

2.4.2 ENERGY BALANCE FOR THE INTERMEDIATE GLASS PLATE 


On the intermediate glass plate, as the depth of water mass 
on the plate increases, the fraction of energy incident on the 
plate decreases. The energy flux terms in the intermediate glass 
plate are as shown in Fig. 2.4. 

The intermediate glass plate being at a higher temperature 
than the glass cover looses heat to the cover while it receives 
energy from the absorbing plate. It also absorbs a portion of the 
incoming solar radiation. The above mentioned energy flux terms 
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for the plate can be expressed mathematically as: 

Solar Energy absorbed by the water mass on the plate=T xi x a 

■g o "w 

Solar energy absorbed by intermediate glass plate=(x I )t a 

g o w g 

Evaporative heat flux from the absorbing plate=Q 

e i 

Convective heat flux from the absorbing plate=Q ^ 

Radiative energy flux from the absorbing plate=Q ^ 

Also there is an energy input to the plate in the form of inlet 
water at temperature T. while there is an energy loss duo to the 
outflow of concentrated brine at temperature T^. Besides this, the 
water vapor that condenses at temperature T 2 leaves the plate at 
temperature 

Hence the energy balance equation of the intermediate glass 
plate is given by 

go w go wg2*^w ic ^el ^cl Vl 

= (M -D ) xCp x(T -T )+D.xCp xCT -T )+Q +Q +0 
22 '^W 2 C 1*^W 2 C ^62 C2 r2 

(12) 

where 

^el" ^l’^ 

^cl" ^2^ 

«rr V o - 

HOLLANDS [7] obtained an empirical expression for the heat 
transfer coefficients in free convection of an inclined layer of 
high aspect ratio, heated from below, as follows; 

““wgo' \go^2'\ 

= 1+ 1 .44x1 l-(1708/RA^g^cos*)j'*'xI 1-1708x8 in(l .84>)^ ^/RA^g^cos^l 


+ l(RA ^ cos4>/5830)*' -1] 

wgc 


( 13 ) 


"'Hie- ‘ . 2 ' V*’. 


where 


X = (|X| + X)/2 


is the modified factor for glass cover since the 
air at the water surface may be considered to be saturated while 
that at the condensing surface may be saturated or 
supersaturated 1191. 

For free convection from a vertical plate, the following 

correlations, given by CHURCHILL and CHU122], have been used. 

* 0.67x l + (0.492/Pr)^'^^®]^^® 

(14) 


Pr = Cp^x 


RA 


weg“ ^ X p X p x(T 2 - 


Simi lar ly , 


Nu^ = 


.=h,H,/K =Nu ^ |„. 


. . . (15) 


wgc 


RA = (gH. Vy.a)[(T,-T_) + (P„ i-Fo>‘Po ^ /(2 . 65xP. -P„ 

1 1 ^ a 12 8,1 2 8,2 1 t s 

where modified factor F 2 is given by 119] 

n 97”? 

F2=0.512I(Tj-T2)/Hj1^-^'*" 

2.4.3 ENERGY BALANCE FOR THE ABSORBING PLATE 


1 


)] 


The energy flux terms for the absorber plate are shown in 
Fig. 2.4. The absorbing plate receives energy from the solar 
radiation transmitted through the glass cover, the water mass and 
the intermediate glass plate. The energy balance equation for the 
absorbing plate can be written as follows: 
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In case of the absorbing plate, it hae been assuaed that, 
except for the glass cover, all parts of the enclosure are 
carefully insulated sc that the heat loss fro. the outside surface 
of the absorbing plate to the anblent «ay be neglected. In 
addition, the absorptivity of water fll. on the absorbing plate 
With blackened wet jute cloths, is assumed to be unity. 


2.4.4 BALANCING OF EVAPORATION RATES 

Considering the second effect, the mass of air transferred per 
unit area per unit time by free convection is given by 


/ Cp X (T - T ) 
az C 2 *^a 2 c 


Therefore . 


T„>«a,' '•'wgc’'VhHeg>''‘we>>«’'r To>'*b 


where 

= 1/(W 
a2 


-W ) 

8,2 C 


=1/{(M„/M^)[P^ “ F P /(P.-F P ) 

W a s,2 t 8,2 C S,C t C B,C 

Combining eq.(14) and (15), we get 
V<V«aV’“^gcV^eg^e>/*b'“’s,s/‘^’t-•’s,^> 


-F P /P.-F P 1 
C S , C t c s ,c 


(17) 


(18) * 


=2.591 (h ^ A +h ^A )/A. JIP /(Pa-P„ ) “F P„ „/(P*-F P )1 

wgc w weg we b s ,2 t s.z cs,c t cs,c 


. . .(19) 


Similarly, for first effect proceeding again as above for D^, we 
get 


D,=2.59 hjlP, ,/(Pt-P, ,) -F.P5,3/«V‘'a*'s,s>‘ ■ 

Thus, there exist five unknown variables Tj, T 2 t T^« Dj end D 2 
which have been determined by solving eqs. (1 1) , (12) , (14) , (19) , (20) 
numerically, as described in the next section. 


30 


2.5 NUMERICAL SOLUTION 

The set of five equations are numerically solved on the 
computer, HP 9000/850 system. The following paragraphs are devoted 
to the explanation of the development of the program [Appendix Al . 

The equations involved are non-linear functions of five 
different variables and the use of a computer facilitates the 
determination of the variables. The program uses the services of 
NAG LIBRARY for the solution of these equations. It uses the 
SUBROUTINE C05NBF(froffl NAG LIBRARY), the purpose of which is to 
find a zero of a system of N non-linear functions in N variables 
by a modification of POVELL’s HYBRID method. It uses a 
Jacobian, calculated by forward difference approximation. It 
requires an external subroutine ,to be supplied by the user, which 
calculates the value of the functions of the N variables 
[Appendix Al . 

Since, hourly and monthly data of temperature and solar 
insolation of NEW DELHI are available from the 1 iterature[2} , 
these are fed as the input data to the computer program. The 
design (width and length of still, chamber height of the two 
stages, weir height), climatic (solar insolation, ambient 
temperature and ambient velocity) and operational (total pressure, 
mass flow rate in the lower basin) parameters are given as inputs 
to the program. 

The initial part of the program is devoted to the calculation 
of various heat transfer coefficients using the initial input data 
and the initial approximation. Also, the various tilt angles of 
the still during different parts of the year are used to calculate 
the various terms occurring in the equations These heat transfer 
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coefficients are calculated by using the values of RAYLEIGH 
NUMBERS yielded by various empirical correlations mentioned 


earl ier . 

These 

heat transfer coefficients 

are then 

used 

to 

calculate 

the 

various heat transfer terms 

occurring 

in 

the 


equations. These equations in the form of non-linear functions are 
then used by subroutine C05NBF to calculate a converging solution 
to the five equations. 



CHAPTER-S 

RESULTS AND CONCLUSIONS 


3.1 INPUT DATA 

To compute the distillate production and the temperatures of 
the various plates of the double-effect still, monthly and hourly 
meteorological data for New Delhi were utilized. This data have 
been utilized in Fig. 3.1 which shows the monthly variation of 
global solar radiation for New Delhi. It can be seen that maximum 
global solar radiation is received in the month of May, where 
after it drops rapidly except for a small rise in the month of 
September. 

3.2 DISCUSSION OF RESULTS 

The temperatures of the cover, the intermediate glass plate 
and the absorber plate and the net distillate production have been 
calculated using the computer program as discussed earlier. 

Fig. 3.2 gives the variation of temperature with the month of 
the year for the cover, the intermediate glass plate and the 
absorber plate. It is seen that the general trend for the 
variation of temperatures follows closely the variation of global 
solar radiation with the month of the yearlFig. 3.1]. Thus, the 
maximum temper, ,»*es are recorded in May. As expected, the absorber 
plate shows the highest temperature and the cover the lowest 
throughout the year. For example, in the month of March, 
temperature calculations for the cover, the intermediate plate and 
the absorber plate are 324 K, 333 K and 345 K, respectively. 



Fig. 3. 1 Average daily values of global solar radiation 
for New Delhi 
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Fig. 3.3 Bhows the variation of distillate outputs for the 

first and second effects separately with the month of the year and 

Fig. 3.4 depicts the total distillate output of the double-effect 
still in comparison with the single-effect still. It can be seen 
from Fig. 3.3 that the output from the second stage of the still 
is higher than that from the first stage. This is due to the fact 

that in addition to the solar radiation absorbed by the 

intermediate glass plate and the water mass on the plate, the 
vapor evaporating from the first stage and condensing on lower 
side of the intermediate glass plate, evaporates an equivalent 
amount of water from the second stage by giving up its latent heat 
of condensation. 

The efficiency of the double effect solar still is defined as 

n = X X (D + D )/I 
12 0 

where D and D are the distillate outputs of the first and second 

t 2 

effects . 

Fig. 3.5 compares the efficiencies of the single and the 
double effect stills. It can be readily seen that the 
efficiency of the double effect still is about 35 to 40 percent 
higher than the single effect solar still. It is also observed 
that the efficiencies of both types of still are higher for higher 
levels of insolation However, the efficiency varies within a 
narrow range of 53.5 to 64 percent for a double— effect still as 

against a range of about 35 to 48 percent for a single effect 

still. It is to be noted that, although, the range of monthly 

variation of solar insolation is quite large, the corresponding 
range of variation in efficiency of both types of solar still is 
small. This is due to the fact that at higher levels of 
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Fig .3 . 4'Coinpar i son of distillate outputs of the single and double-effect stills 
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insolation, the losses are also significant. 

Fig. 3.6a displays the hourly variation of global solar 
radiation on a typical day of March for New Delhi. Fig. 3.6b gives 
the variation of the temperatures for the cover, the intermediate 
glass plate and the absorber plate with the time of the day. Fig. 
3.6c shows the total distillate output for the double effect still 
based on hourly solar radiation data. Figs. 3.6a to 3.6c all 
comprise smooth curves and thus highlight the nearly continuous 
variation of different variables over a period of a day in March 
for New Delhi . 

Fig. 3.7 demonstrates the variation of distillate output with 
ambient air velocity. For a typical day in the month of March for 
NEW DELHI, the ambient air velocity varies between 0 to 10 
m/slll. It can be easily seen that as the velocity of the ambient 
air increases, the net distilled water production rate drops.lt is 
maximum for the still air (i.e. velocity of 0 m/s) and gradually 
decreases as the velocity increases. This is apparently due to 
the fact that, although, the heat losses from th cover increase, 
due to the increased air velocity, resulting in lower cover 
temperatures, there is a fall in the temperature of the 
intermediate glass plate and the absorber plate with the net 
result that the effective difference between the plate 
temperatures decreases , thereby causing a drop in overall 
distillate output. This result is in agreement with the theory 
propounded by LOF et al . , that there should be a drop in the 
distillate production rate with an increase in ambient air 
velocity. The same type of results are reported by YEH I18J in the 
analysis of a single effect still. 


081*8 
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F ig .3 . Sa'Hour I y variation of global solar radiation 
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Fig. 3.8 plots the effect of inlet water temperature on the 
distillate output for a double-effect solar still. There is a net 
increase in the net distillate output due to the increased 
temperature of the inlet water to the lower stage of the solar 
still. This is quite obvious from the fact that the amount of 
energy required for the evaporation of the saline water decreases 
as the inlet saline water temperature increases. 

Fig 3.9 shows the effect of reduced operating pressure on the 
distillate output of the double-effect still. It can be seen that 
as the operating pressure is reduced the net distillate increases. 
It is apparently due to the fact that at a lower operating 
pressure the water evaporates at a much lower temperature thereby 
resulting in higher water mass being evaporated. 

Fig. 3.10 shows the variation of distillate output with mass 
flow rate of saline water in the lower stage of the double-effect 
solar still. The distil Isite output falls as the mass flow rate is 
increased. This is due to the fact that the increased flow rate in 
the lower stage would take away some of the incoming heat, thereby 
resulting in the lowering of the absorber plate temperature and 
the subsequent decrease of the distillate output. 

3.3 CONCLUSIONS 

From the results of the present study, the following 
inferences may be drawn: 

a. A double-effect solar still is about 35 to 40 percent more 
efficient than a single-effect still. 
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Fig . 3 . 8 Effect of inlet water temperature on distillate output 



Effect of operating pressure on disti 1 late output 
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b. The net distillate output for a double-effect still varies froa 
0.403 Kg/m^ hr to 0.886 Kg/m^ hr for various months of the year. 

c. The still has to be operated at the optimal mass flow rate 
corresponding to given still design conditions to achieve higher 
levels of efficiency. This mass flow rate should be such as to 
almost equal the distillate output for the particular month. 

d. Efficiency of the still increases with rise in inlet water 
temperature. This fact helps us to use the hot waste water from 
various thermal plants for the still. It is quite convenient to 
have an inlet water temperature of 303 K for the improved 
performance of the still. 

3.4 SUGGESTIONS FOR FUTURE WORK 

The present study is basically a theoretical investigation 
into the performance of a double-effect solar still. The next 
logical step towards the extension of the present work would be to 
construct and study an experimental double-effect still to verify 
the theoretical results obtained in this study under actual 


conditions . 
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IFAIL = 0 

TOL = SQRT(XO£AJF( ) ) 

TOL = Id-3 

call COSNBF(fcn,n,T,fvec,TOL,UA,LWA,IFAIL) 
wri 16 ( ££*. 9 99 ) 1 , T( 1 ) , T ( £ ) , T ( 3 ) , AD t *36 0 0 0 , AD2*36 0 0 0 
. S70*(AD1+AD£)/AI0 

write (£0 . * ) 1 ,T ( 1 ) , T ( £ ) , T ( 3 ) , AD 1 *360 0 0 , AD2*36 0 0 0 , ( ADI 4-AD£ ) *360 00 

end do 

close(ur{it=21) 
end do 

close (unit =££ ) 
closet unit=£0) 


£ t 0 

FORh/vTt' Temperature Calculations for Month ' , I£, //5X , ' T1 * 


999 


F7.3 
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1 ^ ^ 3,^ Tc *= ',F7.3,/5x, 'Dt • SF5.3,8X,' D2 - ',F5.3, 

2 ,hS 3/) 
end 

subroui ine fcn(n,T,fvec,IFLAG) 


c THIS IS ]:C i->rEFNAL SUBROUTINE WHICH NEEDS TO BE PROVI.DED TO THE« 
c SUBROUTII/* CO.. NET FROfl THE NAG LIBRARV THIS SUBROUTINE CONTAINS THE^ 
c VARIOUS C"UaT]ONS JN IHE FORM OF NON-LINEAR FUNCTIONS AND CALCULATES* 
c THEIR VAi. < ^ lO bt FED TO COSNBF TO OBTAIN A CONVERGING SOLUTION TO* 
c THE G1V£L» FSUmTIONS. * 

^ - . 

iiTiplici: double precision (A - 2) 

c o ‘r> r •:> r- ' c o n - t 1 /' A D . A s . a t , A T a . P T , T i , W > L , p , v , H , A I 0 , Ap h i 

c o fTs :T ■:> r ' •: c ri 1 . 1 r? ^ a p ri i , A D l , A D cl 

c o fi. ? n . c r -- i T ' C F £: , C F c . n 1 . h a , h 2 1 » R A 1 , R A 2 s , R A 2 1 . h c s 

c o A- r. o c O ’ 1 t '• 1 F A I L 

1 n t e . 1 ' . < ^ M I I . n 

d 1 m e F I “T ; Ti ! ^ 3 F . f e c ri ) 

J I = T U ) 

i ( 2 ) 

T c = T ( 3 ) 

2 f(T2 LI Tc) TS = Tc 

1 f ( T1 1. • 12 ) T 1 = T£ 

call 1 c 1 r. a, r d ( T 1 ; I 2 , T c f 


AD I =^2 n1 f 1 f PC ( 1 1 ) / { PT-fps (T1 ) ) ) 

I fp5( 7 2)/(PT-CF2*fps(T2) ) ) 

AD2 ~2 2 9+ ( ^ he t I + h2£ 'i‘ As ) /Ab )*•( f ps ( T2 )/( PT-fps ( T2 ) ) -CFc* 

1 f 1 ! c ) / ( P 1 ^ C? c f p 5 ( T c ) ) ) 

cell c D 1 parf f T 1 , T 2 , T c , qe 2 . qc 2 , qr 2 , qca , qra , qel ; qc 1 , qr 1 > t w, Lw ) 

f V€‘ c ( 1 ) = 0 . 1 2 ^ A I 0 qe2 + qc2 + qr£’-qca*-qra 

fvec<2) = ^0.Frf'AI0)'^Lw+(0.102*AIO)*tw+(3.Od“5)*(Ti*“Tc)+qe1 + 

1 qcl +-qr1- ( (3 0 d-5 ) -AD2 ) * ( T2-Tc ) -AD 1 * ( T2-Tc > “qe2-qc2“qr2 

f uec ( 3 > - 0 7225*A I 0* t W-+ ( 2 . 5d-5 ) * ( Ti - Tc )-qe 1 - qcl - qr t 
1 - < i 2 5d-5 ) -.aD 1 ) *• ( T 1 -Tc ) 

i f < I F I L . e q ^ ^ T h r n 

t e c ( I ) = 0 0 
f ‘ e c ( 2 > 0.0 

1 * ■ € c » 3 " 0 0 

end: f 

r e t u r r, 

end 

subrouTjnr lpcrf(Tl,T?. rc,qe2,qc2,qr2,qca/qra,qet,qct,qr1|tw,Lw) 

c * 

c THIS subroutine CALCULATES THE TOTAL HEAT FLUX CONTRIBUTION DUE TO* 
C CONVECTION rRA^DiATION AND EVAPORATION SEPARATELY. THIS UTILISES THE* 
c HEAT TRANSFER COEFFICIENT VALUES CALCULATED FROM SUBROUTINE CALPARD* 
c AND STILL DESIGN PARAMETERS FROM INPUT FILE. * 


implicit double precision (A - 2) 

common /'con 3 t 1 /Ab/As,At,ATa,PT/Ti,U^L,p,v,H> AI0,Aphi 

common /const2/Asphi/AD1,AD2 ^ 

common / con s t .3 / CF2 , CFc / h 1 , ha / h2T , R A 1 ^ RA2s , R A2t / h2s 

qe2=L( T2)*AD2 

qc2=(h2T*At+h2s*As)*(T2-Tc)/Ab 
qr2=0 . 9*n . 36d-1£)*( ( T2 ) **4- ( Tc ) **4 ) 
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qca^ha* ( Tc^Al e ) 

qra=0 . 95» ( 1 . 36d- 1 £ ) * ( ( Tc ) ♦ ♦4- ( ATa ) *#4 ) 

qetrLdt )*AD1 

qc1=^hU(T1-“T£) 

qr1 = 0 . 96*( 1 . 36d-1£)*( ('T1 )**4-(T£)**4') 
tw=(0. 732-0. 056 t{(p+0.0 084) /p)*log(p+0. 0084) 

1 -( 0 0084/p )+( 1 og( 0 . 0084)-t ) ) 

L = 1 - 1 w 
return 
end 


subrout ine initialise(Acosphi,Asinphi ) 

^ .... 

c THIS SUPPCdTir.? CALCULATES THE VARIOUS SINE AND COSIt4E VALUES OF THE* 
c STILL INCf. JL'.'dION MNGIE AND ALSO THE AREAS OF STILL AND WATER SURFACE* 
c IT USES iHi SUCH AS S-TILL 'INCLINATION ANGLE AND STILL DESIGN* 

c PARAMElERi ♦ 

c ■ ♦ 


1 m p J i ■: I t d u b 1 e p r e c i b i e r» (A “ Z ) 

comrric r / i ons x i / Ab . As , A i , aT a , PT , T i , W , L , p , v ^ H , A I 0 , Aph i 

c o n iT. c r. /'■ c o n s t cv A s p h i , A D 1 , A D 0 

c o /n :n o I . ' c o n s T 3 / C T ? , C F c , h 1 , h , h £ t , R A t , R A £ s , R A £' t , h 0 s 

pi =• c £ / 7 

A p n i = A p h I 1 p 1 / 1 8 0 . 0 
A c o 5 rd) 1 - d c o s { A p h i ) 

ABphx - d r- i n ( A p h 1 ) 

Asinphi-CdsinM .S + Aphi))**1 .6 
A b » U ♦ 1. 

As=*(Ab-^/-sphi )/Aco5phi 
At=Ab/Acosphi 

r e T u r n 
end 

subroutine readver 


c * 

c THIS SUPPQUTINE rEAOS THE VALUES OF INPUT PARAMETERS FROM THE INPUT FILE* 
c AND SUPPLIES THEM FOR CALCULATIONS IN VARIOUS OTHER SUBROUTINES. * 

c ★ 


irr.plicit double precision (A - 2) 

c o Hi n 0 n /const t / A b , A s , A t , A T a , P T , T i , W , L / p , v , H , A I 0 , Aph i 
c o tn n o p /' 1 0 n £. t S / A s p hx ^ ADI , A D £ 

cc'fnmori /cofist3 /CFS, CFC;h1 ,ha,h£t ,RA1 ,RA£s,RA2t ,h£5 
read(cl.'f) U,L^H,p 
r e a d { c‘ 1 . * ) P T , V 

r e t u ’ ' n 
end 


double preci?>ion function fps(T) 

imp* licit double precision (A ~ 2) 

rps = 739 4'^-18 .717*T + 3. 1 595*T*T-'£ . 3 0£5d“£* ( T*T*T ) 
t t8.9St9d-*-4*(T*T*T*T)*7.5d-‘3 

T- e t u r n 
end 


double precision function fplus(x) 

implicit double preci.e. ion (A - Z) 

if (x.GT.0.0) then 

f pi US" X 

else 

fplus=0 

endi f ■ f 

return 

end 

c 

double p. recision function Ka(T} 

implicit double precision <A Z) 

Ka=(2. 4l7£+7 .580d-3*T)d-£ 



end 


ret urn 


c 


c 


c 


double pre<ir*jon function R(T) 

impJi'-it double precision (A •“ Z) 

P - { 1 c'9d3"^ . S14d-3»T+1 . 0583d-5*T^T) 
return 


e no 


douh-ie precision function Cad) 

implicit double precision (A - Z) 

Ca=i 1 . 003S4'»'b.?6d-5*T)d3 
return 


e p o 


dC'Uble precision function nufT) 

implicit double p>recision (A - Z) 

nu* 0 . o - 3^ ( 0£*01 + 0 . 74582^1-0 . 00057 1 7 d { T’t ) 

1 ( 8 9 9^-'P.d~7 ) + (Td^T )-( 6 . 2524d-1 1 

ret u r r, 

e n c! 

c 

d 0 u b i V j> 7' e c 1 e 1 { ! u f u /i c t i o n L ( 7 ) 

implicit double pr'ecision (A - Z) 

L^Ci'SOO £'“£ £975t7-1 . 31 05d-“3*T*T)d3 
ret urn 

end 

^ 

c THIS SUPP''UTIbE CALCULATES THE VALUES OF VARIOUS HEAT TRANSFER COEFFICIENTS 
c AND RAYLEIGH NUMBERS WHICH ARE REQUIRED FOR VARIOUS HEAT TRANSFER TERMS 
c IN THE FUNCTIONS IN FCN 

^ 

subroutine caipard(T1,T2dc) 
implicit, double precision (A “ Z) 

CO mm or /con s t 1 / Ab , A s . At , AT a , PT d i ^ W/ L ^ p , V ^ H , AI 0 ^ Aphi 
common / c o n s t £’ / A s p h i , A 0 1 , A Dd 

common / cor>s t 3/ CFS: . CFc , h 1 ^ ha , hEt , RA1 , RAEs , RAEt , hSs 
common / c c n s t 4 / A c o = p h i > A s i np h i 

CF£=^0 51c't( ( (T1-TE)/H)t.tO.E73) 

CFc=^0 6 06 r( < (72-Tc)/H)*’fO . 119) 

RAEt = vS t >> (H’^HiH)’^RdE)*C(T£)/(nu(TE)<^Ka(TE)*(TE + Tc) ) )♦ 

1 f p 5 ( T£ )-CFc*fps(Tc ) )/( E . 65*PT-f ps ( T£ ) ) >* CTE-273 . 0) ) 

R A 1 - ( £ * 98 ] ^ < H-f H ) ( T1 ) *C( T t )/( nu( T1 ) ^Ka ( T1 ) ^ ( T1 + TE ) ) )♦ 

1 f d 1 ~T£ > - ( ( fps d I )-CF£*fp5 (72) ) / ( £ . 65*PT-“f ps ( T 1 ) ) ) * ( T1 -E73. 0) ) 

RA£s = £^ S‘8t tAco5pna-i^(p>‘P>f‘p)’^R(T£)*C(TE)*(TE*“Tc) 

1 / l tmj ( T 1* i K a ( T £ } t ( T £ + 7 c ) ) 

ha^ I £ . 39 d-5 ) M 5 . 7 + 3 . 8^ v) 

has-- ( Kn d£ > / p ) t ( 0 68 + ( 0 . 67 *RA£b + *( 1 / 4 ) ) )/ 

1 ( 1 f- . ?03" M9/t 6) )t-*'(4/9} 

hi •- i a * 7 1 ) / H ) + ( 1 + 1 . 4 4 + fplus( 1 - ( 1 7 08/ (RA1 *Acosphi ) ) ) 

1 ^ i 1-l7 0 8*^Asinphi/<RA1*Acosphi}) 

E +fplusl (RA1+Acosphi/5830)’«‘*(1/3)*'l )) 

h£t = ( Ka ( r£ ) /H 3* ( 1 + 1 . 44 *fplus( 1-*( 1 708/ ( RAEt ♦Acosph i ) ) )♦ 

1 (1-170StA5inphi/(RA£t >*'Acosphi 3 ) 

E +fp}u£( (RA£t’*‘Acosphi/5830)*^(1/3)-1 33 
rex urTi 
end 


c 
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APPENDIX - B 

CURVE FITTING AND METHOD OF LEAST SQUARES 

Let (Xj, y^) ; i = l,2,...,n be a given set ©f n pairs of 
values X being independent variable and Y the dependent variable. 
The general problem in curve fitting is to find, if possible, an 
analytic expression of the form y = f(x), for the functional 
relationship suggested by the given data. 


If y = a + a X + a X^+, 


rth 


.+ a X’ 
k 


is the K degree polynomial of best fit to the set of points 

(x. ,y . ) i = l,2,..,n the constants a , a , a ..... a are to be 
11 o 1 2 k 

obtained so that 


n 


k.2 




is minimum. 


Thus the normal equations for estimating a^, a^ 




2 k 

are obtained on equating to zero the partial derivatives of E 
w.r.t. a^, a^ , a^ , . , . , a^ separately i.e. 


J) y . s na^ + a^J x. + + a^^ x‘ 


*i ^i ' % *i *i ‘‘k 2) *i 


k+1 


summation extending over i from 1 to n. These are (K + 1) 

equations in (K 1) unknowns a , a , a , a and can be 

O A Z iC 

solved algebraically. 


.2k 
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THERMOPHYSICAL PROPERTIES OF WATER AND AIR 

To analyze the effects of various parameters on the water 
production from a solar still, the model involves the use of 
accurate relations for the thermophysical properties of air and 
water. Using the Least-Squares method described previously the 
following relations for the thermophysical properties of air and 
water have been found to fit the data available from the 
1 iteraturel 1 1 in the range 10 < T°C < 100 with a maximum error of 
0.05 percent. 

a. Density of air 

p„ = 1.2923 -4.514 ♦ 10'® x T + 1.0583 x 10“® x T^ 

O. 

b. Thermal conductivity of air 

K = (2.4172 + 7.580 x lO"^ x T) x lo”^ 

cl 

c. Dynamic viscosity of air 

u = (1.7205 +4.564 x lO"^ x T) x lO”® 

Gl 

d. Specific heat capacity of air 

Cp„ = (1.00324 + 6.76 x 10~® x T) x 10^ 
cl 

e. Specific heat capacity of water 

Cp„ = (4.197 - 8.699 x 10““* x T + 3.686 x lO"® x T^ ) x 10® 

f. Saturated vapor pressure of water vapor 

Py= 739. 4+18. 717xT + 3.1595xT® - 2.3025xl0“® xT® +8.9819x lO'^x T^ 

g. Latent heat of vaporization 

= 2500.2- 2.2975 x T - 1.3105 x lO"® x T® ) x 10® 
w 



